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Abstract 

This article focuses on techniques for acoustic noise reduction, signal filters and source reconstruction. For noise 
reduction, bandpass filters and cross correlations are found to be efficient and fast ways to improve the signal 
to noise ratio and identify a possible neutrino-induced acoustic signal. The reconstruction of the position of an 
acoustic point source in the sea is performed by using small-volume clusters of hydrophones (wlm^) for direction 
reconstruction by a beamforming algorithm. The directional information from a number of such clusters allows for 
position reconstruction. The algorithms for data filtering, direction and position reconstruction are explained and 
demonstrated using simulated data. 
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1. Introduction 

According to Greisen, Zatsepin and Kuzmin 
|1I2| the flux of cosmic-ray protons observed on 
Earth should drop significantly at proton energies 
exceeding some 10 eV due to proton interac- 
tions with the cosmic microwave background. This 
drop in the flux is also known as the GZK-cutoff . 
Recent results from the Pierre-Auger observatory 
support these predictions 13 . In addition to the 
direct measurement of the proton flux, the detec- 
tion of secondary particles resulting from the GZK 
mechanism offers a complementary way to inves- 
tigate this effect. Along with other particles, the 
GZK mechanism produces neutrinos via pion and 
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neutron decays. Due to the low flux of expected 
GZK neutrinos, of the order of 100 km^^year^^ 
[4], huge detectors have to be built. Conventional 
neutrino detectors measure the Cherenkov light 
emitted by secondary particles produced in a neu- 
trino interaction. Such an optical detector would 
comprise an extremely large number of sensors 
because the attenuation length (-^ 60 m (5|) of 
light in water or ice is small. Currently alternative 
methods of detecting ultra-high energy neutrinos 
are being investigated |6]. One of these methods 
is acoustic particle detection, for which an upper 
limit is derived in (7|. Intensive studies are cur- 
rently performed at various places to explore the 
potential of the acoustic detection technique. In 
a recent survey [iS], an overview of these exper- 
imental activities is given. The thermo-acoustic 
model, first discussed by G.A. Askarian in 1957 
II9I10I , predicts the creation of sound waves pro- 
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duced by an instantaneous energy deposition of a 
particle cascade in a medium. The energy deposi- 
tion results in a local rise of temperature causing 
the medium to expand (or contract, depending on 
the medium's properties) which in turn creates a 
sound wave, see Fig. [T] 
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Fig. 2. Simulation of a neutrino induced acoustic pulse for a 
cascade energy of 100 EeV at a distance of 1 km from the 
source. The contributions of the energy deposits at distances 
from the shower axis up to 1 cm and 2 cm and for the whole 
shower are shown. Adapted from 1141 . 

2. Background and signal evaluation in the 
ocean 



Fig. 1. Neutrino-induced sound wave. The typical shower 
length is about 10 m, resulting in a disc-shaped volume, in 
which the particle-induced signal can be detected. 



Test measurements verifying the model for wa- 
ter were performed using proton and laser beams 
|11I12[ and electron beams ITSl . The signal ampli- 
tude is found to depend approximately linearly 
on the temperature and vanishes at the maximum 
density thus demonstrating the thermo-acoustic 
origin of the signal. An exemplary signal shape ex- 
pected for a neutrino-induced hadronic shower is 
presented in Fig. |2l 

The signal is bipolar, with a maximum of the 
power spectrum at about 15 kHz. The attenua- 
tion length of an acoustic signal in water is of the 
order of 1000 m at 20 kHz |!4], yielding a large ad- 
vantage over the optical method for the required 
large detector volumes. 

In this article we describe techniques for data 
processing and source position reconstruction for 
acoustic particle detection. 



The challenge of acoustic particle detection is 
not only to develop sensors sensitive enough 
to detect neutrino-induced acoustic signals. One 
essential problem is to separate the signal from 
the background from a wide range of sources, 
amongst them pressure changes due to current- 
induced turbulences, anthropogenic sounds (e.g. 
ships), wind and precipitation, biological sources 
(e.g. dolphins) and thermal noise. Figure |3] shows 
the power spectrum of the surface-induced noise 
measured in the deep sea compared to the signal 
expectation. While the lower frequencies are dom- 
inated by wind and waves, the higher frequencies 
are dominated by the thermal noise |[T5l|. 

In the interesting frequency band for neutrino 
detection (about 5 — 50 kHz), the wind-induced 
noise decreases with rising frequency, reaching a 
local minimum where the thermal noise starts to 
dominate. The maximum in the frequency spec- 
trum of a neutrino signal is close to the minimum 
of the noise spectrum. Integration of the frequency 
spectrum in the relevant frequency band results in 
a noise level of 10 — 100 mPa. The expected sig- 
nal amplitude is around the same value, for sig- 
nals of a lO^^eV shower at distances of about 100 
m. Therefore methods for noise reduction and sig- 
nal filters have to be developed to enhance the 
particle detection sensitivity; examples are given 
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will have local maxima at times t where the signal 
template matches the recorded data indicating the 
existence of a signal at a time t. Another effective 
method to further improve the signal to noise ratio 
is stacking. Applied e.g. in geophysics for decades, 
this method uses arrays of sensors to identify a 
coherent signal. In the case of imcorrelated noise 
r and a fully correlated signal s in sensors /,/ the 
mean instantaneous signal power is 



with (s,'Sj) 
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while the mean noise reduces to 



Fig. 3. Power spectra from surface agitation and thermal noise 
for different wind speeds, and tfie expected spectrum of the 
neutrino signal (arbitrary units) [15i . 
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with (r,-ry) = ^r^^^y. 
Besides increasing the detection sensi- thus resulting in a signal to noise ratio of 



tivity one also has to understand the signal con- 
tamination resulting from point sources other than 
neutrinos. Of special interest are sources emitting 
neutrino-like, bipolar signals. A prerequisite to ul- 
timately calculate the sensitivity of an acoustic par- 
ticle detector is thus the knowledge of the density 
of sources producing bipolar pulses, and the rate 
of such background signals. 

3. Data reduction 

Since neutrino-induced acoustic signals are pri- 
marily expected in a frequency band of 5 — 50 kHz 
the application of a bandpass filter significantly 
improves the signal to noise ratio. Simulations 
of neutrino-induced cascades with energies of 
^cascade ^ lO^^eV show that a band pass filter 
alone is sufficient to suppress the background to 
a level enabling the trigger of a possible signal, 
e.g. by a simple threshold trigger. Sensitivity for 
lower energies can be achieved by using the cross 
correlatiorl ^ I: 

A{t) = {[r + s](E)s){t) = 

r+T/2 , , , (1) 

dT [r + s]{t + T )s{t ) 

t/2 

where s(t) is the expected signal template and 
(r + s) (t) is the measured pressure amplitude con- 
taining background r and possibly a signal s. A{t) 
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where the sums run over the individual sig- 
nals from all N hydrophones and ( ) represents 
the mean. The -\/N-dependence is the best-case 
scenario for completely uncorrelated noise, but 
clearly indicates the benefit of using a larger num- 
ber of sensors operated in clusters. In geophysics 
this method can be applied passively by reading 
out the sensors in parallel since the duration of 
the signal is large compared to the difference in 
arrival time at the different sensors. In our case 
an active procedure known as beamforming has to 
be used. This method will be discussed in detail 
in section HI In Figs. S] to |6] we demonstrate the 
application of the introduced methods for signal 
to noise ratio improvement. The data for the noise 
and the signal were simulated assuming a sam- 
pling frequency of 200 kHz. The noise spectrum 
was modeled using the data from Fig. |3] for 13 
knotJ^ of wind. A cascade energy of lO^^eV at a 
distance of 200 m and six hydrophones have been 
assumed. For this example the signal arrival time 
is the same for all hydrophones. Figure |4] shows 
the output of a single hydrophone with a signal 
included at about 0.0019 s. Identifying the signal 
is impossible in this case. Applying a cross cor- 
relation significantly improves the signal to noise 
ratio of a single hydrophone as can be seen in Fig. 



^ A technique very similar to matched filtering, described in 

e.g. m- 
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|5l Stacking the output of all six hydrophones im- 
proves the signal to noise ratio to a point where 
a threshold trigger can be applied, see Fig. |6l Ev- 
idently a combination of all introduced methods 
is useful in order to successfully identify a signal. 
The advantage of using clusters of hydrophones 
in order to identify signals is a promising ap- 
proach for future detectors. Furthermore, a cluster 
greatly simplifies the reconstruction of sources, 
which will be discussed in the next section. 
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Fig. 4. Output of a single hydrophone with a noise level 
corresponding to 13 knots of wind and a signal of a 10^' 
eV shower at 200m distance, occurring at about 0.0019 s (see 
marker). 
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Fig. 5. After application of a cross correlation the lower fre- 
quencies are removed. The signal to noise ratio is improved. 
The marker indicates the signal position. 



4. Direction reconstruction 

Reconstructing the position of a point source by 
a small cluster of hydrophones (w 1 m'') becomes 
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Fig. 6. Stacked output of six hydrophones. The signal at 0.0019 
s becomes clearly visible. 

increasingly difficult the further the point source 
is away from the cluster. However, direction re- 
construction is well possible and therefore yields 
the prime observable for a small cluster. It is per- 
formed by creating a sotind intensity plot scanning 
all directions in space (47t). Given an array of N 
hydrophones with coordinates r„ {n = 1,2, ...,N), 
the signals p„ of every single hydrophone will be 
shifted in time corresponding to the difference in 
path length of the sound wave to reach the respec- 
tive hydrophone. Hence every direction in space 
corresponds to a set of time differences Af„ in the 
data. For a direction k, the overall signal at the 
time t is given by 



b(k,t) 



N 

J2 WnVn{i 



(3) 



n=l 



where the Wn represent weighting factors for 
the hydrophones. These factors can be adjusted 
to match the directional sensitivity of the hy- 
drophones. In the following calculations Wn = 1 
was used. The time differences Af„ are computed 
assuming a plane wave [H The algorithm scans 
4/1 with a predefined step size by applying the 
calculated time differences to the data, assuming 
a constant speed of sound. Once all directions 
are scanned, the maximum value of the produced 
output indicates the direction of the signal. Fig- 
ure [7] shows a sample output of the beamf orming 
algorithm. 



^ For point sources this is not entirely true, but a sufficiently 
good approximation if the distances are large compared to 
the dimensions of the hydrophone antenna. 
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Fig. 7. Beam-forming output for a set of simulated data in- 
cluding noise and a signal. The plot shows a clear peak for 
the direction of the signal. The antenna configuration is the 
same as in Fig. |8] 

5. Influence of antenna configurations 

The time delays computed for the beamform- 
ing procedure depend on the relative positions 
of the hydrophones in an antenna. Modifying the 
geometrical arrangement of the hydrophones in 
an antenna will alter its directional sensitivity. In 
this study different geometrical setups of six hy- 
drophones were used (an example is shown in Fig. 

mi. 

Hydrophone Positions 



z[m] 




Fig. 8. The hydrophone positions of one local cluster. 

These geometries were evaluated with respect to 
their directional sensitivity by simulating arrival 
times of a wave front from all directions in steps 
of 5°. Noise samples were generated and signals 
added to the noise corresponding to the time de- 
lays for a given direction. The maximum ampli- 



tude Jmax of the beamforming output divided by 
the RMS of the output was computed for all given 
directions. Figure |9] shows a rather homogeneous 
reconstruction quality for all directions. 




Fig. 9. Maximum pressure amplitude from the beamforming 
algorithm divided by the RMS of the data, for all directions. 

The low /max /RMS values for = 0° and = 
180° are a result of the symmetry of the detector in 
the given coordinate system: If a signal comes di- 
rectly from the top or the bottom of the antenna, all 
(p values for the beamforming algorithm are going 
to contribute with the same value. Increasing or 
decreasing the antenna size by a factor of two cre- 
ates a negligible difference. As long as highly sym- 
metric or two-dimensional geometries are avoided 
every setup will produce similar results. 

6. Source reconstruction 

A simple technique for source position recon- 
struction using a given number of clusters is 
presented. After the successful application of the 
beamforming algorithm each antenna i with a 
detected signal will point to a direction fc,- identi- 
fying where an event came from. Since the cluster 
positions a,- are known, one obtains a line pointing 
to the event for each cluster 

di = ai + njki , e M. (4) 

The source position of a signal is in principle given 
by the intersection point of the lines d,-. Due to un- 
certainties in the beamforming the reconstructed 
source location is where the lines are closest to each 
other. Localising the point of closest approach for 
all lines is realized by calculating the square dis- 
tance from a point s to the reconstructed lines 
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Lf(s) = (s-(fli + «,-fci))^ 
with n,- = (s — ai)ki , 

and minimising the sum of these distances. This 
method was tested for a setup consisting of six an- 
tennas and was found to work well for distances 
up to 1 km from the origin. The test was con- 
ducted as follows: First a source coordinate was 
defined, from where a spherical sound wave was 
emitted, propagating through the water at a con- 
stant speed of 1545 ms^^. As sketched in Fig. [TJ a 
particle-induced soimd pulse does not propagate 
as a spherical wave. However, since the technique 
discussed is independent of the radiation pattern 
a spherical wave was chosen as the most general 
case. A depth dependent velocity gradient was ne- 
glected because of the short distances evaluated 
with respect to the bending radius of about 90 
km ||4]. The arrival times of the spherical wave at 
each single hydrophone were computed and the 
beamforming algorithm used to determine the di- 
rections of the wave onto the storey. Those direc- 
tions were then used to reconstruct the source lo- 
cation. See Fig.[lO]for the geometrical arrangement 
of acoustic sensors and the definition of the coor- 
dinate system. 

Typical errors were less than 10 m for sources 
up to 500 m away from the origin. FigurelTTlshows 
the dependence of the resulting error, i.e. the dis- 
tance of the generated to the reconstructed source. 
For each of the plots in Fig. [TT] two of the source 
coordinates were kept constant while modifying 
the third. The top plot in Fig.[TT]for sources placed 
along the z-axis. Up to 600 m above the sea bed 
the error does not exceed 3 m. For z-coordinates 
above 600 m the error increases due to the angu- 
lar resolution (less than 0.5°) of the beamforming 
algorithm. The plot in the middle is for sources 
placed along the y-axis. The symmetry for the pos- 
itive and negative y-values is due to the detector 
symmetry. The error again increases with distance 
as a result of the angular resolution. The bottom 
plot (sources placed along the x-axis) shows some 
asymmetry resulting from the fact that the clusters 
on the line on the left are higher than the clusters 
on the right resulting in best reconstruction val- 
ues aroimd the position of the line on the left. All 
plots demonstrate that using local clusters for di- 
rection reconstruction and using that information 
to reconstruct a source is a good approach for point 
source location. 
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Fig. 10. The geometry of the simulated acoustic setup. The 
lines are 225 m apart, the vertical position of the stories are 
180 m , 195 m and 305 m for the line on the right and 390 m, 
405 m and 420 m for the line on the left. 



7. Summary and outlook 

This article focuses on techniques for acoustic 
noise reduction, signal filtering and source posi- 
tion reconstruction. For noise reduction, bandpass 
filters and cross correlations prove to be efficient 
and fast means to improve the signal to noise ratio 
and to identify a possible neutrino-induced acous- 
tic signal. The signal to noise ratio can be further 
improved by a factor of \/n {n = number of hy- 
drophones) by the use of local clusters and the 
application of stacking algorithms. Source recon- 
struction using directional information from lo- 
calised clusters is a promising approach. Recon- 
struction errors for sources 250 m from the center 
of the detector are of the order of a few meters, 
increasing for greater distances. The angular reso- 
lution is less than 0.5°. 
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Fig. 11. The position reconstruction error for the cluster setup 
presented in Fig. [lO] The error is the distance from the 
assumed source to the reconstructed source position. Top: 
Sources placed along the z-axis of the coordinate system. Mid- 
dle: Sources placed along the y-axis. Bottom: Sources placed 
along th X-axis. 
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